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ABSTRACT 

In this paper, we consider analytic initial conditions with finite energy, whose complex 
spatial continuation is a superposition of a smooth background flow and a singular field. 
Through explicit calculation in the complex plane, we show that under some assumptions, 
the solution to the 3-D Euler equation ceases to be analytic in the real domain in finite time. 
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1. INTRODUCTION 


One fundamental question in fluid mechanics is whether there exists a velocity field u(x, t) 
and a corresponding vorticity field d>(x,f) satisfying the incompressible 3-D Euler equations 

£>t + u ■ Vdi = w - V u (1) 

V-u = 0 (2) 

V x u — uj (3) 

so that for some time t > 0 these fields are singular when u(x, 0) and tD(x, 0) are initially 
smooth and the energy is finite, i.e. 

J d z x | u(x, 0) | 2 < oo (4) 

The answer to this and related question of enstrophy blow up in the limit of infinite Reynolds 
number have important ramifications in turbulence 1 . 

There have been many investigations in the literature to address this question. Exact so- 
lutions are known 2 for which a smooth flow becomes singular in finite time; the velocity field 
in all such cases do not vanish at infinity and therefore the energy constraint (4) is violated. 
This constraint is important since Childress et al 3 showed that without the energy constraint, 
there are explicit singular solutions of the 2-D Euler equations which is impossible 4,5 ’ 6 for 
flows with finite energy. There have been also been a number of mathematically rigorous 
work'’ 8 on this issue. While there are no definitive theorems one way or the other, a number 
of necessary conditions for finite time singularities have been proved. There have also been 
quite a number of numerical investigations 9-12 as well. Nonetheless, a definitive conclusion 
based on these is difficult to reach. Motivated by physical insight in the analogous mag- 
netohydrodynamics equations, arguments based on an assumed multi-scale expansion were 
presented 13 for the occurrence of finite time singularity for flows with bounded velocity at 
infinity, rather than bounded energy. 

Recently, there have been investigations involving singularities in the complex spatial 
domain. Such studies become relevant to the issue at hand if complex singularities are 
found to impinge the real physical domain later in time. Caflisch 14 simplifies the 3-D Euler 
equation through a so called “Moore approximation”. Within this approximation, a finite 
time singularity is predicted. However, the relevance of this calculation to the actual 3-D 
Euler equation is unclear since the initial data in the real domain is complex and the initial 
value problem is locally ill-posed unlike the actual 3-D Euler equations. Tanveer & Speziale 15 
recently looked at a special class of initial condition that are analytic and real in the real xj, 

^3 domain (x — (xi, X 2 , X 3 )) with a complex spatial plane continuation (i.e. complex xj, 
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x 2 , £3) that includes singularities on a set of points x that are zeros of a a scalar function 
of x. At such points, the analytically extended velocity is finite but vorticity infinite. The 
initial condition is decomposed into a smooth spatial flow and a singular part. Under certain 
assumptions on existence of solutions in this extended complex domain, it was shown that this 
decomposition remains invariant with time, though initial singularities become advected with 
the smooth part of the velocity field (call it the background flow). To the extent that these 
assumptions hold for any class of initial conditions, the results posed no additional restriction 
on the choice of the background flow other than that it is spatially analytic in a neighborhood 
of complex points where the remaining part of the velocity and vorticity fields are singular. 
Indeed for steady background flows, it was seen that the dynamics of the complex singularities 
resulted in sufficient conditions for instability through analysis of three ordinary differential 
equations. Interestingly enough, these differential equations are identical to those found by 
Lifschitz k Hameiri 16 and Vishik & Friedlander 17 , who investigated instability of 3-D Euler 
flows using a WKB analysis for short wavelength disturbances. The latter results have also 
been established rigorously. 

Recently, aside from presenting arguments for finite time singularity for a class of flow 
with bounded velocity but not energy, Bhattacharjee et al 18 applied the complex singularity 
approach of Tanveer k Speziale 15 to a background flow that is spatially smooth in the 
complex plane, but temporally singular. They presented arguments suggesting that if the 
the assumptions made by Tanveer k Speziale 15 are valid, complex singularities can impinge 
the real x domain in finite time resulting in an additional collapsing length scale besides that 
of the background flow. However, no information on the dynamic evolution of the singular 
part of velocity and vorticity is obtained. Further, it is not clear if an initially analytic 
divergence free velocity field in the real domain satisfying the energy constraint (4) can be 
decomposed in the complex x domain into a spatial y smooth background flow that has 
explicit solution and a singular part. 

In this paper, largely motivated by the recent work of Bhattacharjee et al 18 , we present 
analytic initial conditions u(x , 0) in the unbounded real domain satisfying the energy con- 
straint (4) for which the analyticity strip width of solution to (l)-(3) is found to shrink to 
zero in a finite time under some clearly stated assumptions. Since the assumptions above 
have not been verified, the results do not provide a 'dear answer to the question of finite 
time singularity for finite energy smooth initial conditions; nonetheless, we suspect that the 
results will be of value in settling this issue since the assumptions can be checked numerically 
through appropriately designed calculations. 
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2. INITIAL CONDITIONS 


Consider initial velocity field given by 


u(x, 0) = u s (x, 0) + f{d[x, 0]) <?(x,0) , 

where 

- n x t 1 P P 1 , 

U s (X, 0) = ( — -Xj - - x 2 , -X! - -x 2 , X 3 ) , 

m = <?i\ 

d(x,0) = g{ r 2 ) , 

where x = (xi,x 2 ,x 3 ), r 2 = x 2 + x\ + x 2 and 

/ \ _ a 2 (1 + g) 

9 1 + 2s + a 2 (s 2 + s 3 ) ’ 

for some small real and positive constant a 2 and 


(5) 

(6) 

(7) 

( 8 ) 

(9) 


q(x, 0) = [ h(r 2 ) - r 2 ] u,(x, 0) + 2 h'(r 2 ) x x W (x) , (10) 


where 


W(x) 


Hr 2 ) = 

1 

~2 X2X3 


2 g(^ 2 ) 

^(r 2 ) 

1 


+ r" 



( x \ + ^2) 5 


(11) 

(12) 


The initial condition (5) corresponds to an axi-symmetric flow with swirl. Nonetheless, we 
do not make much use of this symmetry since we closely follow the formulation of Tanveer h 
Speziale 15 that was developed for the general 3-D Euler equation. Now, we list below some 
claims that are shown to be true. 


Claim i. The function q[x , 0) is real and analytic everywhere in the real (xi, x 2 , x 3 ) domain. 
It is also analytic in the finite complex x plane except at points where g'(r 2 ) = 0. These 
points are far away from the real domain for small a 2 . 

First notice from (9) that 


g'(r 2 ) 


[1 + 2r 2 + q 2 (r 4 + r 6 )] 2 


1 + 2 a 2 r 2 + 4a 2 r 4 + 2a 2 r 6 j 


(13) 


Notice that the above is nonzero for real x for which r 2 > 0. When a 2 is small, a zero in 
the complex (xi, x 2 , x 3 ) plane occurs only at points far away from the real domain, i.e. for 
I Im Xj| large, atleast for some j. Further, while there can be poles of g(r 2 ), h(r 2 ) defined in 
(11) is clearly analytic at such points. Thus, the function h(r 2 ) is analytic except at points 
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where g'{r 2 ) = 0. The same is clearly true for h'(r 2 ). This establishes claim (i) since 

« s (i, 0) and W(x) are obviously analytic functions of x, i.e. of complex variables Xi, x 2 and 
x 3 . 

Claim ii. The initial condition u{x, 0) is analytic everywhere in the real (x u x 2 , x 3 ) domain 
with singularities in the complex plane determined by the relation g(r ) = 0, g(r ) — oo 
and g'(r 2 ) = 0. 

We note that d(x, 0) = g{r 2 ) is analytic and nonzero everywhere except at points in the 
complex plane when the denominator in (9) vanishes. For small a 2 , zero of the denominator 
in (9) occurs near r 2 = — while the other zeros occur far away from the real domain. 
d(x, 0) is also zero in the complex (x x , x 2 , x 3 ) plane at points where r 2 = - 1 . Raising d 
to a fractional power, we obtain singularities of f(d[x, 0]) at r 2 = —1 (where g{r ) 0) and 

at other complex points where g{r 2 ) = oo. From the decomposition in (5) and claim (i) 
above, claim (ii) follows. 

Claim Hi. The energy constraint (4) is satisfied. 

To show that constraint (4) is satisfied, it suffices to show that for real x with r -* oo, 
it(x, 0) = 0(r -3 ). From the expression of g(r 2 ) in (9) and g'(r 2 ) in (13), it is not difficult 
to show from (11) that for large r, h(r 2 ) = 0(r~ 2 ), h'\r 2 ) = 0{r~ 4 ). Further, from (12), 
W = 0(r 2 ) and therefore from (5), it follows that u(x,0) = 0(r' 3 ). So claim (iii) is 
established. 

Claim iv. The initial condition u(x,0) given by (1) satisfies the incompressibility condition. 
Indeed, each of u,(x,0) and f{d) q{x, 0) is divergence free. 

The divergence condition can be shown to be valid directly by calculating the divergence 
of the right hand side of (5), using the expressions for /(d[x, 0] and (f(x,0) as shown above. 

3. DYNAMICS 

Given that the initial data is analytic in the real domain, it would appear reasonable to 
assume the following, which is listed as 

Assumption (i) 

The initial value problem (l)-(3) with analytic initial condition (5) has an analytic solution 
in the real domain for t in some interval (0, to) for some to > 0. 

Note: If an analytic solution to (l)-(3) does not exist beyond < 0 , then we will have shown 
that the analyticity width shrinks to zero in finite time. In that case, there is nothing more 
to demonstrate. So, we will assume otherwise and let t„ = oo. 
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We now define u s (x,t), hence defined to be the background flow, to be 


u s (x,t) = 


1 -t 


1/3/3 1 

— 2' Tl “ 9 X2 ’ 2 Xl ~ 2 X2 ' X3 


( 14 ) 


corresponding to which the vorticity u s = V x u s is 


^ S ^) 


(3 


1 - t 


[ 0 , 0 , 1 ] 


(15) 


Notice that (14) and (15) satisfy the Euler equations (l)-(3), and are singular at t = 1, 
though there are no spatial singularities; however, since it is unbounded at oo, the energy 
constraint (4) is not satisfied by u s . 

We define d(x, t) to be the solution to 


dt 4“ u s • V d — 0 

satisfying initial condition (8). From the standard method of characterestics, 
d{x{t),t) = d(x{0), 0) = 5(^(0) + ®a(0) + x|(0) ) 
along trajectories determined from 


Solving (18), 
where 


dx _ 
~dt = Us 


r(<) = [r 1 (t),a- 2 (t),x 3 (<)] 


x 3 (t) = 


^ 3 ( 0 ) 
1 - t 


Xi(t) = Xi(0)(l — t) 1 / 2 cos 

^ /n(l - t) 

— x 2 (0)(l — t) 1 / 2 sin 

7} M 1 - t ) 

x 2 (£) = r 2 (0)(l — t) 1 / 2 cos 

-to 

1 

if 

1 1 

+ ri(0)(l — t) 1 ? 2 sin 

^ Ml “0 


(16) 

(17) 

(18) 

(19) 

( 20 ) 

(21) 

( 22 ) 


From (17)-(22), it follows that 

d(x, t) = 9 (j^Tt ^ + x '] + (1 - i) 2 *3]) (23) 

We now decompose velocity u(x,t) and the vorticity u(x, t) = V x u(x, t) into 

u(x, t) = u s (x,t) + f {d[x, t]) q(x, t) (24) 
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c j(x,t) = u a {x, t) + f'(d[x,t]) p(x,t) (25) 

where initial condition (5) is satisfied. Then, in order for the vorticity u(x, t ) and the velocity 
u(x,t) to satisfy (l)-(3), it is clear on substitution that p{x,t) and q(x,t) must satisfy 

P f" f 

Pt + u s • Vp — p ■ — V • q + f(q • V) p + -~(q ■ V) u; s = q(uj s ■ V d) 

J J 

+ (p • V) u s + JJ (u;» ■ V) q + g/(Vd) • (V x q) + fp ■ Vq (26) 

j;V-q + q-Vd = 0 (27) 

xq + Vdxq = p (28) 

The existence of analytic p (x, t ) and q (x, t) in the real domain, satisfying (26)-(28) for t 
in the interval (0, 1) follows from assumption (i) and the decomposition (24) and (25). 

We now introduce another assumption: 

Assumption (ii) 

p(x, t) and q(x,t ) are analytically continuable upto the set of points x 0 (f ), where d(x 0 {t),t) = 0 
Further, each of V q , and Vp is bounded at d = 0. 

Remark: Since each of f and jj are continuous functions of d at d = 0, equations 

(26)-(28) can be seen to have continuous coefficients a; d — 0. So, p(x,t) is expected to 
be continuously differentiable at the set of points x 0 (t), atleast for early time. Moreover, if 
p(x 0 (<), t) is known apriori to be just a locally analytic function of t, then the results below 
will hold for p(x 0 (<), t) for all t in the interval (0, 1), regardless of the assumed boundedness 
in assumption (ii). 

We now define 

A(t) = Vd(x 0 (t),t) (29) 

B(t) = q(x 0 (t),t) (30) 

C(t) = p(x 0 (t),t) (31) 

We denote 

x o (0) = (£,*7,0 

Since d(x o (0), 0) = g(£ 2 + q 2 + ( 2 ) = 0, it follows that 

e + v 2 + C 2 - -i 


From (23) and (29), 


A{ 0 ) = 2g'(-\) [i, 7 , C] 


(32) 
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Further from (10) and (30), it follows that 

B( 0) = 2 A'(-l) f - | v ((’ + -, 2 ) - if C 2 , j f (f 2 + 1, 2 ) - i, C 2 , i C (f 2 + , 2 ) 


It is clear that since d(xo(t),t) = 0, each of /, // /' are zero at d = 0, then with assumption 
(ii), (27) and (28) imply 

AB = 0 (34) 

A x B = C (35) 


It follows that 


C x A = (A- A)B 


Also, from (32), (33) and (35), it follows that 

C(0) = 4 <A— 1) V(— l) (-£( c(f 2 + , 2 )-i,c, j- ,c(f 2 + , 2 ) + iff. 


-f eie + n 2 ) - j« 2 +-? 2 ) 

Furthermore, with assumption (ii), it follows from (26) that 

^ = (dj s -A)B + TC 

where T is a second rank tensor whose elements are defined by 


_ U U c . 

T* = 


In our case, the only nonzero elements of T are 

J ' ^ __ rp rp f rp rp 

In ~~ 2 iT~T)~ l22 ' T 33 ~ r^i ’ ri2 = ~ r21 


l 

2(1 — t) 


From (38), we get 


where S is a second rank tensor whose elements are given by 


Sjk = Tjk + tjklAi 


A ■ A 


given that tjki is the usual Levi-Civita tensor and the Einstein summation convention on 
repeated indices has been used. 
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For t > 0, it. follows from (23) and (29) that 

2<7 # (-l) 


A{t) 


1 - t 


Xl 0 (t), x 2 0 {t), (1 - tf *3o(0 ] 


(43) 


where f 0 (<) = [ x i 0 (t), x 2 0 (t), Z3 0 (0] is determined from (20)-(22) with x 0 (0) - (£, r}, ()• 
Using (20) and (21), it follows 

4/(-l) 


A A = 


1 


[e+v 2 + a -o 3 c 2 


(44) 


Using (43) and (44) in the expression for S in (42), we find that the components of equation 
(41) satisfy the following system of ordinary differential equations: 


dCi 

dt 

dC 2 

dt 


= T\\ C\ 4- T \ 2 C 2 A 


C/5 


= T 2 i Ci + T 22 C 2 + 


dC, 


— T33 C3 + 


C 2 + r/ 2 + (l-0 3 C 2 
C/5 

e + v 2 + (1 - *) 3 c 2 

c/5 


[(1 - tfx 3.(0 C 2 - x 2p (() C 3 \ (45) 

[x 2 „(0 C 3 - (1 - tf *3,(«) Cl] (46) 
[x 5o (t) Cl - x io (t) C 2 \ 


(47) 


dt ' e+T) 2 + (l-f) 3 <2 

Exact solutions to this is possible when ( = 0 in which case each of the components C\, C 2 
is identically zero initally and remains so for t > 0. Using /(—l) = — — *) = 

it follows from (37) and (47) that the third component of C is given by 

3a 2 jd 


C 3 (t) = 


(2a 2 - l) 2 (1 -t) 


(48) 


Therefore for ( = 0, 

c f0 = - 

1 (2a 2 - l) 2 (1 -t) 

It follows from (37), (43) and (44) that in such cases 


0 , 0 , 1 ) 


(49) 


= why* “*■ 01 (50) 

It is clear from (49) and (50) and the definitions (29)-(31) that each of p(x 0 (t)J) and 
q[x 0 (t),t) is singular as t — ► 1~, at the same time as complex singular points 


x 0 (t) = (xi 0 (t), x 2o (t), 0), 

determined from (20) and (21) approaches (0, 0, 0) as (1 - f) 1/2 . 

Since the total vorticity is given by (25), it follows that sufficiently close to these complex 
singular points x 0 (t), 

c o(x,t) ~ YZ~ t I°» °> !) + \ [^(*) • (* - ^o(t))]“ 1/2 C(t) (51) 
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where A(t ) is determined from (43) (with x 3o (f) = 0 in this case) and C(t) is determined 
from (49). It is clear from (50) that the singularities of c5 s and f'(d) p do not cancel out as 
t — + 1~ at So (t). 

In the general case, when ( ^ 0, (45)-(47) do not appear to be exactly solvable. 

Nonetheless, asymptotic considerations suggest that in such cases, all components of C(t) 
generally become singular as t — > 1“. 

4. IMPLICATIONS OF THE CALCULATIONS 

The results (51 ) in the complex plane together with the result that certain class of initially 
complex singularities Xo(0 impinge (0,0,0) in finite time, would appear to suggest that the 
initial value problem (l)-(3) with initial condition (5) becomes singular in the real domain 
in finite time. This would be true if one can assume the following: 

Assumption (Hi) 

If analytic solution to (l)-(3) with initial condition (5) exist for t in (0, 1), then as t — ► 1", 
one or more singularities Xo(t), whose nonzero components are computed through (21) and 
(22), is located on the same Riemann sheet as the real physical domain. 

Note: If assumption (iii) is correct, then the point (0,0,0) at which £ 0 (t) impinges will 
actually be part of the physical real domain and not of another Riemann sheet. 

In the event that assumption (iii) is valid, then the calculations in section 3 show that 
that there must be a collapsing physical length scale near the origin, proportional to (1 — f) 1 / 2 , 
corresponding to the complex singularity distance from the real domain. This is in addition 
to the length scale associated with u s (x s ,t) in (14) that obviously scales as (1 -t). This was 
suggested by Bhattacharjee et al 18 . 

Another point to note about the superposition (25) is that for a fixed real x different 
from 0, as t — * 1 _ , to need not be singular; temporal singularities in u} s can cancel out that 
°f f'{d ) p , as we suspect it will, though we are unable to demonstrate this convincingly. 

The arguments above are not applicable if the solution to (l)-(3) cease to be analytic 
in the real domain for t < 1, as can happen if complex singularities of p, created in the 
dynamical process, impinges the real domain for t < 1. Note that for small a 2 , there exists 
complex poles of g{r 2 ) and hence of d(£,0) near r 2 = -4; from examination of (26)-(27), it 
is clear that singularities of p are expected to be created at such points at t = 0 + . 

Whatever the case, our calculations above show that if assumptions (i)-(iii) are valid, 
then the solution to (l)-(3) ceases to be analytic in the real domain beyond t = 1. It would 
be interesting to check these numerically through schemes that monitor the analyticity width 
as a function of t. Determination of the physical nature of such singularities, if they occur, 
will also be interesting. 
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5. CONCLUSION 

We have considered an initial condition for the three dimensional Euler equation that 
is analytic in the real x domain and contains finite energy. By making three clearly stated 
assumptions, we find through explicit calculations that the solution must lose analyticity in 
finite time. Further study, both analytical and numerical, is likely to shed some light on 
the validity of the assumptions. If these assumptions are true, numerical calculations is also 
likely to suggest if the solution loses analyticity before t = 1 or if it does so at t = 1 in the 
manner shown explicitly in this paper. 
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